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It has been suggested by Timmermans [Phys. Rev. Lett. 87, 240403 (2001)] that loss of fermions 
in a degenerate system causes strong heating. We address the fundamental limit imposed by this 
loss on the temperature that may be obtained by sympathetic cooling of fermions by bosons. Both 
a quantum Boltzmann equation and a quantum Boltzmann master equation are used to study the 
evolution of the occupation number distribution. It is shown that, in the thermodynamic limit, the 
Fermi gas cools to a minimal temperature ksT/fi oc (7i OS s/7coii) ' 44 , where 7i oss is a constant loss 
rate, 7 co n is the bare fermion-boson collision rate not including the reduction due to Fermi statistics, 
and /j, ~ /cbTf is the chemical potential. It is demonstrated that, beyond the thermodynamic limit, 
the discrete nature of the momentum spectrum of the system can block cooling. The unusual non- 
thermal nature of the number distribution is illustrated from several points of view: the Fermi 
surface is distorted, and in the region of zero momentum the number distribution can descend to 
values significantly less than unity. Our model explicitly depends on a constant evaporation rate, 
the value of which can strongly affect the minimum temperature. 

PACS numbers: 



I. INTRODUCTION 

Evaporative coolinghas proven able to obtain degen- 
erate Fermi systems ]J0|lSEIi0ll]- As polar- 
ized fermions cannot undergo s-wavc collisions, it is nec- 
essary to sympathetically cool with another species or 
spin state. In contrast to the case of bosons it has 
been suggested that Fermi systems are highly sensitive 
to loss [l(| • In this article, we will investigate this ques- 
tion for the case of sympathetic cooling by an ideal zero 
temperature Bose gas, in order to identify the funda- 
mental limit imposed by loss. Sympathetic cooling via 
a degenerate Bose gas is indeed used in several present 
experiments 0, S S ■ The lowest experimentally ob- 
tained temperature to date is T/Tp ~ 0.05 [g, where Tp 
is the Fermi temperature. 

Various theoretical groups have considered improved 
cooling schemes [Til ], including optimizing the evapora- 
tion rate 01 1 using different trapping frequencies 

or using laser rather than evaporative cooling [la . 
[la . Il7j . We will restrict our investigation to a simple 
model of sympathetic cooling of a gas in a box, in which 
it will be shown that a minimum temperature arises nat- 
urally as a result of loss of particles. The discrete na- 
ture of the energy spectrum of the system can also be 
a limiting factor. It will be demonstrated, from several 
points of view, that the occupation number distribution 
is non-thermal in a non-trivial sense. "Temperature" and 
"chemical potential" will therefore be defined based on 
the total number and energy of the fermions, rather than 
on the equilibrium nature, or lack thereof, of their distri- 
bution. 

Our presentation is organized as follows. In Sec.^ we 
physically motivate our idealized model in the context of 
present experiments. In Sec. IIIII a quantum Boltzmann 



equation is derived under the assumption that the den- 
sity operator is gaussian in the fcrmionic field, which per- 
mits the use of Wick's theorem, to study the evolution of 
the mean number distribution and the temperature. In 
Sec. II VI this equation is investigated numerically in a dis- 
crete system and both numerically and analytically in the 
thermodynamic limit. In Sec. the probability distribu- 
tion of the occupation numbers is examined without the 
assumption of Wick's theorem: a quantum Boltzmann 
master equation [l8L Hsj based on the secular approxi- 
mation [23, [21j is derived and simulated via Monte Carlo 
methods. Finally, in Sec. IVll we conclude. 



II. A MODEL FOR SYMPATHETIC COOLING 

In order to understand the limits of sympathetic cool- 
ing imposed by particle loss, we use an idealized theoret- 
ical model. This model contains the following assump- 
tions. Firstly, the bosons form a perfect reservoir: they 
are non-interacting, at zero temperature, are not signif- 
icantly reduced in number during the total observation 
time, and when excited to non-zero momentum states are 
removed from the system by evaporation. Secondly, the 
Fermi gas is non-interacting: the fermions are all in the 
same spin state so that there is no s-wave contribution 
to their interactions, and, in the low temperature regime 
which will be considered, the p-wave contribution is neg- 
ligible. Thirdly, fermion-boson interactions arc modeled 
by a contact potential proportional to the s-wave scatter- 
ing length. Fourthly, the system is enclosed in a three- 
dimensional box with periodic boundary conditions with 
sides of length L Xl L yi L z : the volume shall be denoted 
as L 3 = L x L y L z . Fifthly, the evolution of the system is 
described by discrete evaporation time steps of duration 
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t e , the end points of which represent a full removal of the 
bosonic particles in the modes with non-zero momentum: 
l/t e may be interpreted as the evaporation rate. Sixthly, 
a constant fermion loss rate is introduced. 

In present experiments it is commonly assumed that 
the Fermi and Bose gases are cooled down while remain- 
ing thermalized at a finite and common temperature. 
This requires a sufficiently large collision rate between 
bosons and fermions as well as between bosons, in com- 
parison to the evaporation rate. In our model, each time 
an excited boson is created by interaction with the Fermi 
gas it is removed sufficiently rapidly so that, even if an 
interacting Bose gas was considered, the bosons would 
not have time to thcrmalizc. Consequently, the Fermi 
and Bose gases are never at the same temperature. It is 
in fact advantageous to maintain the Bose gas at a tem- 
perature much smaller than that of the Fermi gas: it is 
only in this case that all collisions between fermions and 
bosons are efficient, in the sense that they decrease the 
energy of the fermions and therefore lead to cooling. 

In this respect, our model is not intended to closely 
represent current experiments; rather, the goal of this 
study is to explore the fundamental cooling limits due to 
loss in an ideal system. However, we note that in an ac- 
tual experiment the regime considered in our model may 
be obtained if thermalization of the bosons is avoided 
by a sufficiently strong evaporation rate. We also note 
that cooling of fermions by a nearly pure condensate was 
realized in a recent experiment || . 

III. THE QUANTUM BOLTZMANN EQUATION 

Consider the Hamiltonian H = Ho + V, where Hq is 
the kinetic energy of the Fermi-Bose mixture and V is 
the interaction energy, such that 

h 2 h 2 , h 2 n 2 

k q 

F =~^ 51 fr q ^qCk'Ck<5k'+q\k+q, (2) 
k',k,q'/q 

g(t) = g sin 2 (tit /t e ), (3) 

where & (b>) refers to the creation operator, subject to 
the usual Fermi (Bose) commutator relations, for a sin- 
gle fermion (boson) of momentum k (q). The coefficient 
go = 2nh 2 a s /m r , a s is the s-wave scattering length for 
fermion-boson interactions, and m r is the reduced mass. 
The choice of a continuous time dependence for the inter- 
action potential in Eq. ( avoids strong non-adiabatic 
effects. In contrast, an abrupt switching on and off of 
the interaction potential is a source of heating, leading 
in particular to divergence of the mean kinetic energy in 
the rate equations to follow. A Feshbach resonance could 
be used to achieve a time-dependent coupling of the form 



suggested by Eq. (H) |H H3 . Note that the case q = q' 
has been discarded in Eq. (0, as it gives a contribution 
gNftN-p/ L 3 which has no effect on the dynamics. 

In order to apply perturbation theory, it is required 
that t e be sufficiently small. For example, in the ther- 
modynamic limit in the Fermi degenerate regime, it is 
necessary that 

teJcollT/Tp < 1 , (4) 

where 

7coll = |crn-B«F , (5) 
o 

cr = Aira 2 is the cross section for scattering between a 
boson and a fermion, n-Q is the bosonic density, and v-p 
is the Fermi velocity, where vf = Ti^im-p) 1 ^ / m-p in this 
spin-polarized system, with ?ip the density of fermions. 
The factor of T/Tf in Eq. (@J is due to Pauli blocking. 
The factor of 3/8 in Eq. (JSJ has been included to account 
for the reduced interaction strength due to the choice 
of the temporal profile of g(t). Note that 7 co ii is time- 
dependent for a non-zero loss rate, since the density of 
fermions decreases. Consider the number operator 

C^4c k . (6) 
The mean value of £ at time (n + l)t e may be written 
(0[(n + l)t e } = Tr[p(nt e )A(nt e )], (7) 

where 

A(i) = 17+ [(n + l)t e *-t]C U[{n + l)t e «- t) , (8) 

and p(nt e ) and U are the density operator and the evo- 
lution operator in the interaction picture. The operator 
A satisfies the equation of motion 

ihj t A=[V(t),A(t)}, (9) 

with the "final" condition A[(n + l)t e ] = C- This is equiv- 
alent to the integral equation 

r(n+i)i e d 

A perturbativc development of A may then be obtained 
by iteration of Eq. (|10|) in powers of V as follows: 

r (n+l)t c dt _ 

A(nt e )^(+ ^[{,V(t)} 

Jnte IK 

One must then calculate the expectation value of A(nt e ) 
with respect to the state of the system after the evapora- 
tion step at time nt e , which is defined in the interaction 
picture by 

p{nt e ) = p F (nt~) ® \N B : q = 0)(A B : q = 0| , (12) 
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p F (nt e ) = Ti B [p( n t e )] 



(13) 



The fact that, after evaporation, all remaining bosons are 
in q = was used. Note that the number of bosons has 
been assumed to remain approximately constant during 
the evaporation process, in keeping with the assumption 
of a perfect reservoir. Depletion of the reservoir can only 
decrease the cooling efficiency. 

The mean occupation number of the single particle 
state with momentum k is defined as 



iV„(k) = Tr[cj c c k p(nt e )] 



(14) 



As shown in App. [XJ the time development of Eq. Hll|) 
plus the use of an approximate Wick factorization leads 
to an approximate rate equation for the occupation num- 
bers which iteratively describes the development of the 
system in evaporation steps of period t e : 

N n+1 (k) = (l- 7loss i e )A„(k) 

+ ^P(k'^k)7V„(k')[l-Ar„(k)] 

k' 

- ^7>(k -> k')[l - N n (k')]N n (k) ,(15) 



where 



7>(k->k') = 



h 2 L^ 



dt g(t) exp(iwi) 



hu> = 





h 2 k' 2 

2mp 



fi 2 (k-k' 
2m B 



h 2 k 2 

2mp 



(16) 



,(17) 



(18) 



and the integer n is the previous number of iterations. 
In the right-hand side of Eq. I|15f) . the second term rep- 
resents the sum of probabilities of moving a fermion into 
state k while the third term is the sum of probabilities 
of moving a fermion out of state k. A loss term with 
a constant rate has also been introduced, under the as- 
sumption 7ioss^e <C 1. This describes, for example, col- 
lisions with background gas particles present in experi- 
ments. Equation l|15|) is the central result of this section 
and the basis of further study in this article. Its validity 
is subject to the necessary condition that the probability 
of departure from mode k after an evaporation cycle is 



k'^k 



k')[l-iV„(k')] « i 



(19) 



for all populated levels k. A similar condition must hold 
for the probability of arrival. 

When t is expressed in units of the evolu- 

tion of N n (k) is ultimately governed by three continu- 
ous dimensionlcss parameters, 7i OS s/7coii, teE-F(0)/K, and 
fint e / '(mpL 2 ), as well as the ratio of masses a = mg/mF, 
which is fixed for a particular experiment. As the goal 



of this work is to study the ultimate limits of sympa- 
thetic cooling, the minimum temperature in units of the 
chemical potential shall be studied as a function of these 
parameters. Temperature and chemical potential are de- 
fined with respect to the total number of fermions and 
the total energy, as given by the standard sums over k of 
N n (k) and (h 2 k 2 /2m F )N n (k) |53 , respectively, not with 
respect to the equilibrium nature, or lack thereof, of the 
number distribution. 

One may ask if there are higher order effects on 
Eq. (|15fl which arc important. The third order term 
of Eq. (|llf) has a vanishing contribution; however, 
the fourth order term contains several physical effects. 
Firstly, it contains a correction to the Born approxima- 
tion for the scattering of a single fermion with a single 
boson; this correction is small provided that k F a <C 1. 
This is equivalent to the weakly interacting regime, since 
fcp oc n^/ 3 . Secondly, a boson may interact with a 
fermion and leave the condensate, undergo a subsequent 
interaction with a second fermion, and enter a final mo- 
mentum state q". There are then two possibilities: if 
q" = 0, bosonic stimulation occurs, and the contribution 
is proportional to [7Jb|(?(uO| 2 ] 2 ; in contrast, the sum over 
all final states q" ^ has a contribution proportional 
to n F n-Q\g{^)\ A ■ The former, which represents effective 
interactions between fermions mediated by the bosonic 
reservoir, has been studied elsewhere [28|, and is here ne- 
glected. The latter is smaller than the former by a factor 
of Np J Nb < 1. 



IV. STUDY OF THE QUANTUM BOLTZMANN 
EQUATION 

In the following, Eq. i|15|) is studied with three different 
approaches. In Sec. lIV A*l the evolution is investigated in 
a discrete system via numerical integration. In Sec. lIV 51 
the thermodynamic limit is taken and a second numerical 
study is made. Finally, in Sec. IIV CI the thermodynamic 
limit is treated approximately and analytically under the 
assumption of an equilibrium Fermi distribution. 



A. Time Evolution for a Finite System 

We first consider the evolution of a discrete, finite sys- 
tem which evolves according to Eq. 1)150. There are two 
distinct regimes of t e . When the half- width 2tt /t e of the 
function g(u>) in Eq. I|17|) is large compared to the spacing 
5E between the values of Tilo in (|15fl . the discrete nature 
of the spectrum of the system does not play a significant 
role: this we term the continuous regime. A calculation 
of a typical SE is presented in App. [Bj Using the result 
of this calculation, the continuous regime may be defined 
explicitly by 



SE t e 
2nh 



8 (EfU\ (a + l) s 
3 V KN F J 4a 2 



< 1 



(20) 
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where a is the mass ratio ms/mf. If, in addition, iV(k) 
varies slowly with respect to the momentum level spac- 
ing, one may furthermore take the thermodynamic limit, 
as shall be elucidated in Sec. HVBl For SE t e /2irh S> 1, 
g{ui) is no longer well-resolved and the discreteness plays 
a strong role: this we term the discrete regime. 

Figure shows a study of Eq. H15|) under varia- 
tion of the two central parameters 7i OS s/7coii(0) and 
t e Ep(0)/h for the case of a 7 Li - 6 Li mixture, as in 
Refs. HQ. Typical experimental values of 7ioss/7coii(0) 
arc 8/3 x 1CT 2 [Fig. [Tfa)] to 8/3 x 1CT 3 [Fig. □»]. 
The value 7 io S s/7coii(0) = 8/3 x 1(T 4 [Fig. Die) Which 
could be reached by using a Feshbach resonance [22, 123] 
to augment the scattering length a s and thus the col- 
lision rate 7 co n cx a 2 , was simulated as well. Simula- 
tions of Np = 10 2 (dashed curve) and Np = 10 3 (long 
dashed curve) fermions were performed in a nearly cubic 
box with incommensurable sides; for computational rea- 
sons, a cube was used for Np = 10 4 (dot-dashed curve). 
One clearly observes the continuous regime to the left 
hand side of each plot, where the temperature is domi- 
nated by the evaporation rate, as explained in Sec. II V 51 
The optimal temperature is achieved in the vicinity of 
t e Ep/Ti ~ 10 2 . Towards the right hand side of each plot, 
the minimum temperature rises due to the blocking of 
cooling by the discrete nature of the spectrum; this ef- 
fect weakens for larger numbers of atoms for a given value 
of Ept e as is apparent in Eq. (|2l)[l . One may observe 
the blocking in Eq. I|17|) . where, for values of the en- 
ergy difference hu> h/t e , g{uj) becomes small and the 
interaction is effectively reduced. Finally, for values of 
t e Ep/Ti > 10 4 perturbation theory is no longer applica- 
ble, as 70011^/^(0) = 8/3 x 10~ 5 was chosen, in keeping 
with typical experimental conditions |24j| . 

Figure [21 shows the effect of t e on the occupation 
number distribution in energy space. In the continuous 
regime the distribution clearly depends on energy alone, 
so that in the thermodynamic limit A^(k) = N(k). The 
distribution is not, however, an equilibrium one, as is 
more easily observed in the thermodynamic limit (see 
below). In contrast, in the discrete regime the distri- 
bution is fully k-dependent. Careful observation of the 
figure shows regular holes in the energy spectrum: this is 
a natural result of the quantization of the box. We note 
that, in all regimes investigated, the essential feature of 
a Fermi surface, though distorted, persists. 



B. Time Evolution in the Thermodynamic Limit 

In the thermodynamic limit, as defined explicitly in 
Sec. IIV Al the sums in Eq. (|15|l may be approximated 
using the standard continuum limit 




E 



(21) 
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FIG. 1: Shown is the minimum value of ksT/fi as a func- 
tion of the time step t e Ep/h for a 7 Li - 6 Li mixture and 
fi7coii/£ F (0) = | x 1(T 5 . The ratio 7 i oss /7coii(0) is 8/3x 1(T 2 , 
8/3 x 10~ 3 , and 8/3 x 1CT 4 in panels (a), (b), and (c), re- 
spectively. The case of a finite system, Eq. (1151 . is shown 
for 10 2 (dashed line), 10 3 (long dashed line), and 10 4 (dot- 
dashed line) fermions, and the thermodynamic limit, Eq. (1231 . 
is shown as a solid line. For the finite system, to the right 
side of the plot, in the discrete regime, the cooling is blocked; 
to the left, in the continuous regime, the minimum temper- 
ature is dominated by the evaporation rate, according to 
Eq. and rises as t e Ep/h decreases. The data obtained in 
the thermodynamic limit are independent of t e except where 
the evaporation-limited temperature is larger than the loss- 
limited temperature. Note that the actual data points are 
represented by open circles [25ll . 
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FIG. 3: Time evolution of the number distribution of a Fermi 
gas sympathetically cooled by a perfect Bose reservoir for 
two experimentally important cases in the thermodynamic 
limit: (a) 7 Li - 6 Li and (b) 23 Na - 6 Li. Shown are the initial 
state of (ksT/ [i)(0) = 0.7 (thin line), the state at which the 
temperature reaches a minimum (thicker line), and a later 
stage where the gas has heated (thickest line). The hole in 
the distribution near k = is particularly evident in panel 
(b), where the loss processes balance evaporative cooling, as 
described by Eq. Here [7k,ss/7cou(0), (kBT/fj,) mill \ = 

[8.27 x 10" 3 , 0.081], and [1.25 x 10~ 2 , 0.089] for panels (a) 
and (b), respectively, and h'Jcoii/ Ef(0) = | x 10~ 5 . 



FIG. 2: Shown is the mean occupation number for JVf = 1000 
6 Li atoms in contact with a 7 Li reservoir in a finite sys- 
tem. The initial state (plus symbols) is a thermal distri- 
bution of k B T/fi = 0.324; h-y coll /E F (0) = § x 10~ 5 . Af- 
ter an evolution time (solid circles) such that ksT/fi passes 
through a minimum according to Eq. I|15^ with the parameters 
7io SS /7coii(0) = 8/3 x 10" 3 and (a) t e E F /h = 10 2 (continuous 
regime) (b) t e E F /h = 10 3 , and (c) t e E F /h = 10 4 (discrete 
regime) , it is clear that in the continuous regime the distribu- 
tion, though non-thermal, depends only on |k|. In contrast, 
in the discrete regime the distribution is fully k-dependent. 



As illustrated in Fig. [21 in this limit 

N(k) = N(k) . (22) 

The continuum iterative equation for the time evolution 
may then be written as 

N n+1 (k) = (1 - 7ioss t e )N n {k) + ^ / dk ' k ' 2 

n Jo 

x |jV n (A)[l - N n (k')] J ^ du\g(oj)\ 2 -(k~ fc')} ,(23) 
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where g{u>) and uj are defined as in Eqs. (|17|l and 1181) 
with 



k • k' = kk'u . 



(24) 



Two effects limiting the choice of t e are implicit in 
Eq. (J2HJ- The first, which appears when t e is small, is 
the width of the interaction function g(u>): as k^T ap- 
proaches the width Anh/t e in the course of cooling, the 
sharp decrease in N n (k) from unity to zero, typical of 
a Fermi or Fermi-like distribution, is no longer resolved, 
and the transfer of momentum from fcrmions to bosons 
ceases to have any effect on the temperature. This gives 
an absolute minimum temperature of 



te 



(25) 



The superscript "evap" refers to the fact that l/t e is the 
evaporation rate. The second effect, which occurs in the 
limit of large t e , is the validity of perturbation theory, 
according to Eq. J3J. Since this limitation is imposed by 
our use of perturbation theory it is not fundamental to 
sympathetic cooling [26|. Therefore, in practice, within 
the context of our model, one should choose an evapora- 
tion rate such that 



T~p Tt? E~p 

— < t e E F /h < — ^~ 

T T 7coll h 



(26) 



We have verified that the results are independent of val- 
ues of t e which satisfy Eq. I|26fl. as illustrated by the 
plateau in Fig. ^ 

It is convenient to begin with No(k) in the form of 
a Fermi distribution (though other initial distributions 
were studied, with the same qualitative results). In the 
following simulations, k^T/ ' [i = 0.7 was taken as a start- 
ing condition. In Fig. [3] is shown the evolution of the 
occupation number distribution resulting from Eq. 1231) 
with a choice of t e satisfying Eq. iffiH) and 7i OS s/7coii(0) = 
8.27 x 10~ 3 and 1.25 x 10~ 2 for the experimentally rel- 
evant cases of 6 Li- 7 Li and 6 Li- 23 Na, as shown in panels 
(a) and (b), respectively. Close inspection of the fig- 
ure shows that the distribution is non-equilibrium: the 
thermal tail is missing, and, as can be seen by making 
a fit to a Fermi distribution (not shown), the rise from 
N n (k) = 0.5 towards unity with decreasing k is less sharp 
than that of a Fermi distribution with the same total en- 
ergy and number of fcrmions. There is also a hole near 
k = 0, which is difficult to see in Fig. |3{a) but appears 
strongly in Fig.[3{b). One may observe the existence of 
this latter feature directly from Eq. Il2.'{[) as follows. 

The evolution of N n (k = 0) may be approximated by 
assuming that N n (k) varies slowly near the origin, which 
allows one to replace N n (k') with iV n (0) in Eq. J53J). The 
condition that the value of N n (0) increase is then 



7io 



7coll 



< 



4C 
3^ 



[l-iV(0)] 



aEpt c 



1/2 



(l + a 2 ) 



1 



|1_q,|3/2 (l +a )3/2 



(27) 



where a = ms/m-p and C = 1.860266 .... It is therefore 
directly apparent that for sufficient loss rates the number 
distribution has a hole at k = 0. Moreover, since a fac- 
tor of [1 — N n (0)] enters into the condition, the number 
distribution never reaches unity at k = and the dis- 
tribution is never a Fermi distribution. It was observed 
numerically that this feature extends up to a finite fc, 
the width of which varies dynamically and increases as 
m^/m-p takes a value largely different from unity. The 
deepest point occurs at k = 0, and the maximum in time 
of N n (0) is given by replacing the less than or about equal 
to sign with an equal sign in Eq. I|27|) . The evolution of 
the hole for a 23 Na- 6 Li mixture is particularly apparent, 
since a ~ 23/6 is far from unity, which reduces the value 
of the right hand side of Eq. (|27|) for a fixed value of 
N n (0), as illustrated in Fig.^b). 

The effect of the choice of 7i OS s/7coii on the minimal 
temperature is shown in Fig. The power laws 



h T\ ( \ 0.436 

= 0.659 [-^] +5.6xl0- 4 (28) 



/ mm 

knT 



7coll 



0.445 

= 0.621 ( — 1 +7.34xl0~ 8 (29) 



/' 



7coll 



for 7 Na- 6 Li and 23 Na- 6 Li, respectively, were found over 
the range 8/3 xl0~ 5 < 7i oss /7coii < 8/3xl0" 2 . Note that 
the constant offsets in the above are negligible over the 
fit domain. Although the distributions shown in Fig. |3| 
are not Fermi distributions, the step-like feature makes 
a temperature and chemical potential, as defined by the 
total number of particles and the total energy, a mean- 
ingful measure of the shape of N n {k). Moreover, as iVp ot 
and E F ot are weighted by k 2 and fc 4 respectively, the hole 
near k = has little effect on them. 

Finally, Fig. [S] illustrates the time evolution of k^T '/ '/i 
for the parameters 7i OS s/7coii(0) = 8.27 x 10~ 3 and 
t e Ep(0)Ti = 10 3 . The figure is divided into three time 
regions which show three phases of the evolution: cool- 
ing, the achievment of a minimal k^T/fi, and heating. 
These three regions were observed in all simulations for 
which 7ioss/7coii *C 1, and were independent of t e Ep/h 
within the constraints of Eq. (|26fl . 



C. Analytical Prediction of the Degeneracy for a 
Fermi Distribution Ansatz 



The following concerns the thermodynamic limit alone. 
In the limit in which 



h/t e < k B T , 
one may use the approximation 



(30) 



(31) 



where S(uj) is a Dirac delta distribution. Substituting 
Eq. (|31|) into Eq. (|23|) . and defining time continuously 
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FIG. 4: The numerically determined maximum Fermi degen- 
eracy as a function of loss rate for 7 Li - 6 Li (open circles) and 
23 Na - 6 Li (black diamonds) mixtures in the thermodynamic 
limit, according to the quantum Boltzmann equation. One 
finds (fcBT/V)min x (7ios B /7coii(0)) ' 44 . The dashed line shows 
the analytic prediction of (feT//j) m i„ cx (7ioss/7coii(0)) 1/3 
based on the assumption of a Fermi distribution. The dif- 
ference is due to the non-thermal nature of the quasi-static 
mean occupation number distribution, an example of which 
is illustrated in Fig. EH Here /2/y co ii/£f(0) = § X 1CT 5 . 



according to 



N(k,t) 



Nn+l - Nn 
t. 



(32) 



the iterative rate equation reduces to a first order partial 
integro-differential equation with an integration over k' 
alone: 

*(M) = -WV(M) + 1=# 

o 2irn 
dk' —N(k',t)[l - N(k,t)] 




dk' — [1- N{k',t)]N(k,t) 
k k 



(33) 



where a = ms/m-p. Here the functional form of the 
limits of integration was determined by integration of 
the delta distribution over the solid angle. Note that in 
the limit in which k — > 0, N(k) — — 7i oss iV„(A;), provided 
that ttt-f ^= rriB- In the limit as m B — > tof = Tn, the 
limits of integration of Eq. 1)33(1 simplify to [k, oo] and 
[0, A;], respectively. In this case the 1/k in the integrand 
is not regulated and the expression diverges as k — > 0, 
save in the case where [1 — N(k)] — > 0. 

An analytical model of the time evolution of the tem- 
perature can be developed based on a Fermi distribution 
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FIG. 5: Time evolution of the degeneracy in the ther- 
modynamic limit as determined by numerical evolution of 
Eq. 123H (solid line), and the analytical prediction of Eq. 13711 
(dashed line) based on a Fermi ansatz. The parameters are 
7ioss/7coii(0) = 8.27 x 1CT 3 , (fe B T//x)(0) = 0.7, t e E F (0)h = 
10 3 , and hf coU /E F (0) = § X 10" 5 . 



ansatz, namely, 
N a (k,t) 



CX P{[^T-M(t)] AbT(*)} 



(34) 



For simplicity, the case m = tof = % is considered. 
Two equations for the unknowns T(t) and /i(i) are ob- 
tained by multiplying Eq. I|33|) by k° and k 2 and inte- 
grating over k: 



dt 



JV* ot 



""floss- 1 »F 

L 3 



(27T) = 



h 2 k 2 

2m v ; 



(35) 
(36) 



where iV F ot and -E F ot are the total number of fermions and 
total energy of the fermions, respectively. In the degener- 
ate regime where k B T/ fi <C 1, one may obtain an approx- 
imate evolution of k B T j [i from Eqs. Q35J) and l|36(l . One 
neglects terms of order exp(— fi/k^T) in the right hand 
side of Eq. I)36|l and uses the low temperature expansions 
of Np ot and £ F ot H3|. Details are given in App.O One 
finds 



d 
dt 



k B T 



■ r )7T 



.-Tic 



12C(3) 



Tcoll(*) 



k B T 



(37) 

which clearly shows the separate contributions of heating 
due to losses and cooling due to elastic collisions. The 
fact that the cooling term is proportional to T 3 has a 
simple physical interpretation. Each collisional process 
occurs at a rate JcoII^bT / (J. due to Pauli blocking. It 
involves a fraction k B T / \i of the total number of fermions. 
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The energy transferred to a boson per collisional event is 
of order k B T. Therefore the collisional term in dE/dt is 
proportional to T 3 , from which follows Eq. I|37)l. 

In the limit 7ioss/7coii "C 1, which is in fact the exper- 
imental case, one may distinguish three different stages 
in the evolution of k B T '/ 'fx. In the first stage, cooling 
dominates and it decreases according to the power law 



(*) 



1 



6C(3h- 2 7 coll (0)*+(/i/fc B T)(0) 



(38) 



Note that, in the case where 7i oss = 0, Eq. I|38|) holds 
indefinitely. In the second stage, after a time t m i n , k B T/fj, 
arrives at a minima, given by 



fc B r 



7loss 



1/3 



15C(3)7coii(0)_ 



(39) 
(40) 



Note that Tiossimin "C 1, so that a very small frac- 
tion of the atoms have been lost when k^T/p achieves 
its minimum. In the third stage, heating manifests as 
an adiabatic increase in k B T/fi, obtained by neglecting 
d(k B T / p,) 2 / dt in Eq. (|33(l and thereby replacing 7 co ii(0) 
by 7coli(i) = 7coli(0) cxp(-7i oss t/3) in Eq. figft. The evo- 
lution of k B T j p continues to increase adiabatically up till 
a characteristic time given by 



^-nonadiabatii 



3 . I Tcoll 
In ' 



Tie 



(41) 



Note that, at this evolution time, k B T / ' p is on the order 
of unity, and the above analytical treatment ceases to be 
applicable. 

The number A TVs of bosons lost from the reservoir dur- 
ing the first and second stages of cooling, i.e., up till the 
time at which the minimum k B T/p is achieved, may be 
estimated simply from Eq. (|38|l . One integrates the rate 
of production of excited bosons oc NF(0)(k B T/p) 2 'y co ii(0) 
over time, and finds the approximate relation 



AiV B oc N F {0) 



knT 



(0). 



(42) 



This corresponds to the initial number of fermions active 
in the cooling process. 

In Fig.^is shown the minimum temperature as a func- 
tion of loss rate. The difference between the analytical 
prediction of Eq. l|3^|) (dashed line) and the numerical re- 
sults of Eq. (|23|l (solid line) highlight the non-equilibrium 
nature of the actual mean occupation number distribu- 
tion. However, the qualitative features are the same 
for both the analytical model and the numerical simu- 
lation. Figure shows the evolution for the parame- 
ters 7io SS /7coii(0) = 8.27 x 1CT 3 and (k B T/p)(0) = 0.7. 
The three stages of cooling, achievement of a minimum 
temperature, and heating, are clearly observable. As 
7ioss/7coii(0) <C 1, the time scale of the first stage is small 
compared to the third stage. 



V. BEYOND THE BOLTZMANN 
APPROXIMATION: THE MASTER EQUATION 

The quantum Boltzmann equation is a closed equa- 
tion obtained after use of Wick's theorem to replace the 
mean value of occupation numbers with the product of 
their mean values. Such an approximation applies when 
the probability distribution of the density operator is 
nearly Gaussian [Tsl llfij . and neglects correlations be- 
tween modes. In the following, the exact probability dis- 
tribution of the occupation number shall be treated by 
deriving a master equation for the fermion density opera- 
tor in the limit of weak coupling. Specifically, g must be 
small enough so that the probability to have more than 
one boson excited out of the condensate during t e is neg- 
ligible. In the Fock basis, the density operator is char- 
acterized both by 'populations' and 'coherences', that is, 
by the diagonal and off-diagonal matrix elements, respec- 
tively. Subsequently, the secular approximation will be 
used to derive an equation for the evolution of the dis- 
tribution probability of the occupation numbers. This 
approximation applies in the regime in which the evo- 
lution rate of the populations is much smaller than the 
Bohr frequency of the coherences [2(j, which allows one 
to derive closed equations for the populations. 



A. Derivation of the Quantum Boltzmann Master 
Equation 

Just after a measurement of state of the bosons, but 
before excited bosons have been removed from the system 
via evaporation, the fermion density operator may be 
written 

M(n + l)t e ] = ( n b ■ 0\p[(n + l)t e ]\N B : 0) 
+ - 1 : 0; 1 : q\p[(n + l)t e ] \N B - 1 : 0; 1 : q) (43) 

where a sum over bosonic Fock states has been taken. 
By choosing go sufficiently small one may indeed neglect 
the possibility of exciting more than one boson out of the 
condensate. The first term in Eq. I)43|l represents zero ex- 
cited bosons, and the second term a single excited boson 
of momentum hq. The combined bosonic and fermionic 
density operator p[(n + l)t e ] is given by the action on 
Eq. (|12f) of the evolution operator U(nt e ) from time nt e 
to time (n + l)t e : 

p[(n+l)* e ]= (44) 
U(nt e ) [p F (nt-) ® \N B : q = 0)(N B : q = 0|] W{nt e ) . 

One then expands the evolution operator U(nt e ) to sec- 
ond order in the interaction potential V using standard 
time-dependent perturbation theory, which results in the 
second order expansion of Eq. (|44[l. and thus Eq. (|43fl . 
One obtains the following operators which act on the 
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FIG. 6: Evolution of the total energy of the fermions for 
a finite system for one Monte-Carlo realization of the exact 
occupation number distribution under the secular approxi- 
mation (see text) according to Eq. 1571 . (a) 7i O ss/7coii(0) = 
8/3 x 10~ 4 : the loss events (in this case 6) are marked by 
arrows. Clearly cooling dominates the decrease in the energy 
at early times while loss plays a stronger role at later times, 
(b) 7ioss/7coii(0) = 8/3 x 10~ 2 : The upper curve shows the 
evolution of the energy, and the lower curve shows the evo- 
lution of the number of atoms. At early times, the slope of 
the energy is compounded of the effect of cooling and a trivial 
decrease of the Fermi energy tx exp(— 2-y\ oss t/3) due to atom 
loss. In both (a) and (b), (feT/p)(0) = 0.324, jV F (0) = 100, 
t e E F /h = 10 3 . 

fermions alone: 
Cq{nt e ) = 



A(nt e ) 



(N B -l:0,l:q\V(t')\N B :0), (45) 



(n+l)ie ft' dt n 

t ih 



c io~' 
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FIG. 7: The numerically determined maximum Fermi degen- 
eracy as a function of loss rate for N = 100 (dashed line) and 
N — 1000 (solid line) 6 Li fermions in contact with 7 Li (open 
circles) or 23 Na (black diamonds), according to the quantum 
master Boltzmann equation. The value of t e E F /Ti which gives 
the minimal temperature has been chosen. Arrows indicate 
the long term cooling limits for 7i oss = for 7 Li- 6 Li and 
23 Na- 6 Li mixes, respectively. The minimal temperature is 
much higher in comparison to that predicted by the quan- 
tum Boltzmann equation in the thermodynamic limit, as was 
shown in Fig. 2] (note scale change). 



{N B :0\V(t')V(t")\N B :0) 
C q may be written explicitly as 



C q (nt e ) 



hk 2 



2mp 



Tiq 2 
2niB 2m F 



(46) 

(47) 
(48) 



The explicit form of A is more complicated, but may be 
derived by substituting Eq. J5J into Eq. |@BJ|. C q orig- 
inates from the second term of Eq. I|43|) in which one 
boson is excited. It originates from a first order per- 
turbative expansion of U, but appears in two factors of 
Eq. (|44|l and therefore gives a contribution of order g 2 in 
the evolution of the density operator. A originates from 
the first term in Eq. (43) with the evolution operator 
expanded to second order in go. It contains, in partic- 
ular, an effective interaction between fermions mediated 
by the bosons. This results in the master equation for 
the fermionic density operator 

PF[{n + l)t e ] = pF{nt e ) + A(nt e )p F (nt e ) 
+p F {nt e )A\nt e ) + ^2 Cq{nt e )p F {nt e )Cl(nt e ) . (49) 

We now proceed to apply the secular approxima- 
tion [2^. To illustrate the details, the contribution of 
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Cq is described explicitly. Evaluating the last term in 
Eq. gjjj, 



[C q/ 5 F C q ](nt e ) 



E 1 

ki,k 2 



g(ui )g* (w 2 )cj ti „ q c kl p F (n* e )cJ C2 c k2 _ q , 



(50) 



where cJi^ is defined as in Eq. I|48|) . with k replaced by 
k 12 . The typical evolution rate 7 CV oivc 01 the fermionic 
density operator is proportional to g 2 ,. In contrast, the 
Bohr frequencies u)\ — uji do not depend on go. There- 
fore the oscillating exponential in Eq. (|50|l may be ne- 
glected for sufficiently small g , as its effects averages to 
zero when averaged during l/7 C voivo- An estimate for 
Tevoive may be made based on the assumption of a ther- 
mal distribution in the thermodynamic limit, as defined 
by Eq. ||13J): 



""/evolve 



Tcoll-^F 



TV 



(51) 



where SN-q <C 1 is the mean number of excited bosons 
during one cycle of duration t e . The minimal Bohr fre- 
quencies arc given by 



for which, in Eq. (|50() . one keeps only the terms with 
lui = lu 2 - Assuming the box lengths squared to be in- 
commensurable, one finds 



or kic 



(54) 



for each a € {x,y,z}. The former case, which corre- 
sponds to the excitation of a boson in the plane orthog- 
onal to a, we neglect. The existence of this solution is a 
consequence of the separability of the degrees of motion 
along x, y, z in the box. One could consider an alternate 
model for the box in which this separability is lifted to 
justify its being neglected ^5|. There therefore remains 
the sole condition k x = k 2 . 

Having applied the secular approximation, if p F is ini- 
tially a statistical mixture of Fock states, then it remains 
one for all times. Defining the occupation number prob- 
ability distribution Q„({n k }) by 



p F (nt e ) = Qn({n k })\{n k })({n k }\ , (55) 
{™k} 



min(|wi - u) 2 \) 



1 



hp(E F ) 



Airmkp 



(52) 



where p(E F ) is the density of states at the Fermi sur- 
face. Therefore the condition for validity of the secular 
approximation is 



E F 



T 
TV 



< 1, 



(53) 



where 



|{n k }> = n(4)"Vc 



(56) 



and each n k G {0, 1}, one obtains the equation of motion 
for Q: 



Qn+i({n k }) - Q„({n k }) = E kl #k 2 n°L 6 {-n^ii 1 - "-k 2 )Qn({n-k}) + «k 2 (l - n kl )Q n ({n k + c5 k , kl - 5 k .k 2 })} 

- (Ek 7lo SS tenu ) Q({nk}) + [Xko 7losste(l - «k )] <3({"-k + ftc,k }) , (57) 

I 



with 

^=^1 + h ^ k ^ _ U?i . (5 8 ) 
2mp 2me 2m F 

Here the loss term has been added in by hand under the 
assumption that 7i OS sie-^F <S 1- The condition that the 
number of bosons excited during t e be much smaller than 
unity may be written 

J2 P(ki^k 2 )n kl (l-n k2 )<l, (59) 

ki,k 2 

for typical configurations {n k }, where V is defined as in 
Eq. (|16f) . This is to be contrasted with the much weaker 



condition of Eq. I|19fl obtained in the Quantum Boltz- 
mann equation, where the sum is over only one momen- 
tum. 



B. Monte Carlo Numerical Study 

The continuous time version of Eq. Ij57|) . where Q n +i — 
Q n = t e Q, was studied numerically via Monte Carlo sim- 
ulation. All data for N F (0) = 10 2 and N F (0) = 10 3 in 
Fig. ^ was re-evaluated, with a mean over 100 realiza- 
tions of Q made for each data point. The qualitative 
features of the evolution of the mean number distribu- 
tion were found to remain the same: its non-equilibrium 
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nature, as illustrated for the finite system in Fig. and 
for the thermodynamic limit in Fig. [3J and the three evo- 
lution stages of rapid cooling, achievement of a minimum, 
and slow heating with a quasi-static distribution, as illus- 
trated for the thermodynamic limit in Fig. [5j Two exam- 
ples of a single Monte-Carlo realization are illustrated in 
Fig. |SJ The minimum temperature for these simulations 
is achieved at 7 co ii(0)t ~ 20. Before this time the evolu- 
tion of the energy is dominated by cooling, rather than 
losses: in Fig. HJa) this is directly apparent, as the loss 
of individual fermions is marked by arrows; in Fig.EJb), 
where the loss rate was higher, the slope is seen to be dif- 
ferent from that which results from the mere evolution of 
the total number of fermions, the latter of which is deter- 
mined according to (A F ot ) = N^ ot (0) exp(-j loss t) . After 
this time loss dominates, as is apparent in both panels: 
the greater part of the changes in the energy occur at the 
same time as loss of a fermion. 

However, quantitatively the agreement between the re- 
sults of the quantum Boltzmann equation and the quan- 
tum Boltzmann master equation depends both on the 
number of atoms and on the evaporation rate. For 100 
atoms, the minimum temperature predicted by the mas- 
ter equation, as shown in Fig. [7] is as much as 50% higher 
as for the same parameters as in Fig.^ This deviation is 
strongest for small loss rates, where the cooling efficiency 
is limited by the blocking mechanism due to the discrete 
nature of the spectrum, as discussed in Sec. IIV Al For 
1000 atoms, the agreement is very good for the loss rates 
of Figs.^a) andQIb) while for Fig.^c) there is a 15% 
increase in the temperature. The power laws 



A* / mm 



/ \ 0.503 

0.861 ( — ) + 1.44 x 10~ 2 (60) 

\7coll/ 

\ 0.518 

= 0.855 1 — ) + 0.993 xl0~ 2 (61) 

7coll / 



for 7 Li- 6 Li and 23 Na- 6 Li, respectively, were found over 
the full range of data shown in Fig. [7| for N = 1000. 
Comparing to Eqs. (|29|l and 1|29|1 . which were obtained in 
the thermodynamic limit with the quantum Boltzmann 
equation, one observes that the exponents of the power 
law are similar, whereas the constant offset is substan- 
tially different. Recalling that the data of Fig. [7| resulted 
from an optimization over t e , this suggests that even in 
the absence of loss there is an absolute minimal temper- 
ature of /cbT//i ~ 10~ 2 . This was validated for 1000 
atoms by Monte Carlo simulations with a zero loss rate, 
as indicated by arrows drawn along the left hand y-axis 
of the figure: for 7 Li- 6 Li, mm(k B T/n) = 1.34 x 10~ 2 ; for 
23 Na- 6 Li, mm(k B T/fi) = 0.98 x 10" 2 . 

Therefore, by taking the actual probability distribu- 
tion for the occupation numbers into account, i.e., by not 
assuming near-thermal equilibrium according to Wick's 
theorem, it is found that the minimum temperature in- 
creases when loss does not dominate the cooling. This 
again highlights the non-thermal nature of this system. 
In contrast, when N-p is large and the loss rate is in the 



experimental range of 10~ 2 7 co ii to 10 -3 7 co n, the secu- 
lar approximation shows that Wick's theorem is a valid 
assumption and one may simply use the quantum Boltz- 
mann equation. 



VI. CONCLUSION 

We have investigated the effect of heating caused by 
loss of atoms on the minimum temperature that may 
be achieved in a sympathetically cooled Fermi gas. The 
model of sympathetic cooling that we have used is cycli- 
cal and consists of a sequence of time intervals during 
which fermions are coupled to a zero-temperature ideal 
Bose gas via binary atomic interactions. At the end of 
each time interval the excited bosons are removed from 
the system by evaporation. The length of the time in- 
terval is short enough that the bosons do not come into 
thermal equilibrium with the fermions; this is in contrast 
to present experiments in which the bosons and fermions 
are always in equilibrium with one another. The cooling 
is balanced by a constant loss rate of fermions, which, for 
example, can be caused by collisions with background gas 
in the experimental apparatus, and can lead to heating, 
as shown in Ref. ^j- The combination of cooling and 
heating are modeled at several levels of theoretical ap- 
proximation. 

First, a quantum Boltzmann equation describing the 
evolution of the mean occupation number distribution 
was developed under the assumption that the fermion 
density operator is nearly Gaussian, i.e., that Wick's 
theorem may be applied. The overall minimum temper- 
ature, which was found to be best obtained in the ther- 
modynamic limit, was observed to follow the power law 
(fc B T//u) min ~ 0.65(7i oss /7coii) ' 44 , so that {k B T/^j,) min < 
0.03 for 7ioss/7coii _ 10~ 3 , where /i ~ k B T? is the chem- 
ical potential of the fermions, j co i\ is the bare fermion- 
boson collision rate not including the reduction due to 
Fermi statistics, and 7i oss is the constant fermion loss 
rate. This value of 7i OS s/7coii is easily achievable in 
present experiments, in particular by using a Feshbach 
resonance |23l | . The number distribution was observed to 
have a distorted Fermi surface and a hole near k = 0. 

A second theoretical perspective was developed based 
on the secular approximation to a master equation, with- 
out the use of Wick's theorem. Monte Carlo simula- 
tions of the resulting quantum Boltzmann master equa- 
tion showed that in the limit of experimentally reason- 
able values of 7i OS s/7coii and the number of fermions 
N-p, the first theoretical approach is indeed valid. How- 
ever, for values of 7i OS s/7coii tending towards zero, the 
master equation shows a substantially higher tempera- 
ture. In the most extreme case studied of 100 atoms and 
7ioss/7coii = 3 x 10~ 4 , this increase is 50%. 

A possible extension to this work is to add a harmonic 
trap and/or to include interactions between bosons. The 
assumption of a perfect Bose reservoir is reasonable when 
the speed of sound is much smaller than the Fermi veloc- 
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ity [29j, as is the case for weakly interacting condensates. 
Although a harmonic trap may change the predicted min- 
imal temperature, the qualitative results of this study 
are expected to be correct even for a non-uniform system 
such as is found in present experiments on Fermi-Bose 
mixtures. 
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APPENDIX A: EVOLUTION OF THE 
OCCUPATION NUMBERS 

In the following, we will calculate an approximation 
to the variation of the expectation value of the number 
operator 



C( 



(Al) 



from time nt e to time {n + l)t e , that is, during the time 
interval of duration t e between two successive measure- 
ments of the state of the bosons. 

One begins with the approximate evolution of the num- 
ber operator obtained using second order perturbation 
theory for A, Eq. I|ll|) . One takes the expectation of 
Eq. I|ll|) with respect to the density operator Eq. (|12fl . 
All the bosons at time nt e are in the state with zero 
momentum. In order to take the expectation value first 
with respect to the bosons, it is convenient to rewrite the 
interaction potential V as 



V(t)= ]T WW(*)& q iq 

where the W's are purely fcrmionic operators: 
5o (t) 



(A2) 



Wq'q(t) 



L 3 



e^v.^oA+q^+q' , (A3) 



with 



hk' 2 hq 2 hk 2 hq 2 



■,'2 



^k'q'.kq — 



2mp 2me 2mp 2mj 



(A4) 



Expanding the commutators in (|12|) and using the fact 
that V does not contain terms with q = q', one obtains 
the following matrix elements and their expressions: 

(JV B : 0|y(t)|JV B : 0} =0 (A5) 
(N B : 0\V{t)...V(t')\N B : 0) 

= JV B ^ Wo,q'(*) ■ ■ ■ WV, (*') (A6) 
q'#o 



where (...) may contain an operator acting on the 
fermions alone. These results may be interpreted physi- 
cally as follows. The action of V on a pure Bose-Einstein 
condensate creates a state with A^b — 1 ground state 
bosons and a single excited boson with a non-vanishing 
momentum q', since the terms with q = q' have been ex- 
cluded from the expression of V, as apparent in Eq. J5J). 
The resulting excited state of the bosons is orthogonal 
to the initial state, so that the term of Eq. I|ll|) linear in 

V has a vanishing expectation value. A second action of 

V gives a non-zero contribution to the expectation value 
only if the excited boson is scattered back into the con- 
densate. 

There are also terms where the factors involving the in- 
teraction potential appear in reverse chronological order, 
such as V(t')V(i). These terms are Hermitian conjugates 
of the terms in chronological order so that the final result 
reads 

(Co)[(n + l)t e ]-<Co>(nt e ) = -^E / dt ' 



(n+l)t e 



dt ([Co, Wo, q '(t)]W q /, (t')} + cc. (A7) 



where the expectation value in the right hand side 
is taken with respect to the fermion density operator 
/Of (nt e ). 

Finally, one evaluates the commutator in (|A7|) . The 
identity 



; k £k o' 4' a k] = (Sko.k'C^Ck - S koM cl,c ko , (A8) 

which is a direct consequence of the fermionic anticom- 
mutation relations, is used. Observing that the conserva- 
tion of momentum imposes k' = k + q' in the expression 
for Wo,q', one obtains 



9o(t) 
L 3 



ikoCko-q'e^o -^-^'* 
Cko+q'Ckoe 40 ^' ^'*) ■ (A9) 



[C0 ) Wb,q'W] = 



Multiplying this expression by Wq',0 gives fourth degree 
equation in the fermionic creation/annihilation opera- 
tors: to obtain a closed equation for the occupation num- 
bers one performs a crucial factorization approximation 
based on the Wick contraction rule. This constitutes the 
weak point of the present approach, which was explored 
by a more systematic treatment in Sec. V. As the system 
is spatially homogeneous, the mean value of the product 
of a creation operator and an annihilation operator of 
different momentum states vanishes. One is left with 



(Ck Cko-q'C k _q'Ck) 



&,k JVn(ko)[l-JVn(ko-q / )] 

(^k +q'^koC k _q'Ck) — 
<5k,k +q'[l - N n (k )]N n (k + q') 



(A10) 
(All) 
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where the occupation numbers N n (k) are defined by 
Eq. 114|) and the fact that q' ^ has been used. Observ- 
ing that 

(n+l)t e An+l)te 

dt 1 / dt goi^goit'y^'-^ + c.c. 

nt e J t' 

= I.9M| 2 ,(A12) 

where g{y>) is defined in (|17f) . one obtains the desired 
identity l(T3|) . 



APPENDIX B: LEVEL SPACING IN THE 
QUANTUM BOLTZMANN EQUATION 

Consider a given single particle level of wavevector k 
in the box. In the quantum Boltzmann equation (|15|l . 
this level is coupled to all the other levels k'. We wish 
to estimate the mean distance between the values of the 
corresponding energy mismatches %ui given by Eq. (|18|1 . 
One may then define the density of these values of fiu) by 



Since |g(o;)| 2 is centered in to = 0, one can restrict the 
density of the Tiw's to E = 0. Furthermore, we approxi- 
mate the discrete sum in <r(0) by an integral: 



equilibrium Fermi distribution, respectively Substitut- 
ing Eq. ||22Jl into the right hand side of Eq. , 



4 ot W = -71oss£f 



dx 



X |a;[l - N a (x,t)} J yiV (j/,t) 
- x N a (x,t)£^-[l-N a (y,t)}^ , 



(C5) 



where 
T = 



L_\ (l n B g 2 m B \ ( 2mk B T " 

2tt) \2m) \8 2wh 3 J \ h 2 



(C6) 

the integration variables x = (h k 2 ) / (2mk B T) , y = 
(H 2 k' 2 )/(2mk B T), and N a {x,t) = N a (k,t), N a (y,t) = 
N a (k',t). The last term in Eq. Ip^. 

h = ~-r dxxN a (x,t) dy[l-N a (y,t)}, (C7) 
4 Jo Jo 



(Bl) may be integrated by parts: 



y ,.00 />oo 

J 2 = - - jf [1 — N a (x, t)} J dy yN a (y, 



i). (C8) 



This reduces the integral in Eq. (|C5|) to 



a(0) 



Jd 3 k'5{hu). (B2) 7 =4"/ da: / 



dy[l-^(z,i)]i\r (l/, *)(*-!/). (C9) 



This integral can be calculated exactly by using spherical 
coordinates and integrating first on the polar angle, then 
on the modulus k'\ 



Making the substitutions \ = x ~ ^H^bT), ip = y — 
fi/(k B T) and replacing N a according to Eq. i|34[l . one 
obtains 



L\ rusk 
1 2irl 2h 2 



a-1 
a + l 



(B3) 



where a = itib/itif. Taking the typical value k ~ kp, 
using the relation Np = (L/2ir) 3 4'Kk F /3 and setting 
SE = l/er(0), one obtains Eq. 



APPENDIX C: EVOLUTION OF THE 
TEMPERATURE 

In the following, Eq. i|37[) . which describes the time 
evolution of k B T / fi under the assumptions 



h/t e < k B T , 
m = mp = m B , 
{k B T/fi) 2 « 1 , 
N(k,t) = N a (k,t), 



(CI) 
(C2) 
(C3) 
(C4) 



is derived from Eqs. I|35(l and (|36J) . These four assump- 
tions correspond to the use of Fermi's Golden Rule, equal 
masses of fermions and bosons, high degeneracy, and an 



T 



X l x # ( e ^ + l)(et + i) 



(CIO) 



The key to resolving Eq. (|C10|) is to take ~^/(k B T) — > 
—oo in the lower limit of the first integral. Note that this 
is consistent with assumption (|C3|1 . In this case, 



1 4 
2T 



/ = C(3) = 1.20206 • 



(Cll) 



where £ is the Riemann zeta function [30l | . 

What is the error involved in this approximation? 
Defining the neglected portion of Eq. HC10|) as 



J = 



dX J x ^ + l)(e-x + l) 



(C12) 



it is apparent that the second integral is evaluated over 
€ [x < —/i/(k B T) <C — l,oo]. The leading contribution 
of this integral therefore gives 



d ^WTT) 



I 2 

2 X 



(C13) 



14 



Substituting Eq. HC13|) into Eq. (|CJ12|> . the error is 



J- -- 



fc R T 



Therefore 



4 ot W--7io SS i?r-^TC(3) 



-(9 



(n/k B T) e 



The high degeneracy expansions 

3 



N tot 



L\° . 1 (2mn\* 2 



2- 



4?r 



2\ h 2 J 3 



(C14) 



(C15) 



7T 

1 + T 



(9 



(C16) 



7T! tot 



L\° h z . 1 /2m^\ 5 2 



2tt 



-47T- 

2m 2 







fc R T 



(C17) 



may be easily developed from the treatment of Ref . [27| • 
Sub stitut ing Eq. QCfty into Eq. (gSJ and Eq. (fcT7|) into 
Eq. (|C15|I . one may solve Eq. (|35fl for fi(t) and substitute 
the resulting expression into Eq. (|36J) to obtain the final 
result, (ksT/ fi)(t), as shown in Eq. I|37|) . 
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